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ABSTRACT 

We use data from the ARGO cosmic microwave background (CMB) anisotropy 
■ experiment to constrain cosmogonies. We account for the ARGO beamwidth and 

calibration uncertainties, and marginalize over the offset removed from the data. Our 
derived amplitudes of the CMB anisotropy detected by the ARGO experiment are 
smaller than those derived previously. 
\ We consider open and spatially-flat-A cold dark matter cosmogonies, with 

clustered-mass density parameter fio m the range 0.1—1, baryonic-mass density 
parameter Q# in the range (0. 005-0. 029)/i -2 , and age of the universe to in the range 
(10-20) Gyr. Marginalizing over all parameters but Qq, the ARGO data favors an 
open (spatially-flat-A) model with Qq = 0.23 (0.1). However, these numerical values 
O are model dependent. 

At the 2 a confidence level model normalizations deduced from the ARGO data 
are consistent with those drawn from the UCSB South Pole 1994, MAX 4+5, White 
Dish, and SuZIE data sets. The ARGO open model normalizations are also consistent 
with those deduced from the DMR data. However, for most spatially-flat-A models the 
DMR normalizations are more than 2 a above the ARGO ones. 
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1. Introduction 

Ganga et al. (1997a, hereafter GRGS) developed a general method to account for experimental 
and observational uncertainties, such as those in the beamwidth and the calibration, in likelihood 
analyses of cosmic microwave background (CMB) anisotropy data sets. In conjunction with 
theoretically-predicted CMB anisotropy spectra, this method has been used to account for 
beamwidth and calibration uncertainties in analyses of the Gundersen et al. (1995) UCSB South 
Pole 1994 data, the Church et al. (1997) SuZIE data, the MAX 4+5 data (Tanaka et al. 1996; 
Lim et al. 1996, and references therein), and the Tucker et al. (1993) White Dish data (GRGS; 
Ganga et al. 1997b, 1998; Ratra et al. 1998). 

In this paper we present results from a similar analysis of the ARGO CMB anisotropy data 
set, but now consider a much larger range of cosmological-model parameter space. The 1993 flight 
of the balloon-borne ARGO experiment resulted in two data sets. One was from data taken in 
the direction of Hercules (de Bernardis et al. 1994a, hereafter deB94), and the second consisted 
of data taken in the direction of Aries and Taurus (Masi et al. 1996). In addition to the CMB 
signal, the Aries and Taurus data set is known to have a significant foreground dust signal in the 
shorter wavelength data with milder contamination in the longest wavelength data (Masi et al. 
1995, 1996). To use the Aries and Taurus data to investigate CMB anisotropy one must carefully 
model this foreground dust contamination (Masi et al. 1996). In this paper we therefore only 
consider the deB94 Hercules data. De Bernardis et al. (1993) and Palumbo et al. (1994) describe 
the ARGO telescope and detectors. De Bernardis et al. (1994b) discuss the constraints that follow 
from the ARGO Hercules data on some cold dark matter (CDM) and other cosmological models. 

ARGO data were taken in four wavelength bands centered at 0.5, 0.8, 1.2, and 2.0 mm 
(deB94) . While there is dust contamination in the shorter wavelength data, the longest wavelength 
data is thought to be pure CMB anisotropy (deB94). Following deB94 we use only the 2 mm data 
for our CMB anisotropy analyses here. 

The FWHM of the beam, assumed to be gaussian, is 52', with a 3% one standard deviation 
uncertainty. While observing, the beam was sinusoidally chopped with a half peak-to-peak chop 
amplitude of 0.9° on the sky. 63 points were observed with the payload performing azimuth 
step and integrate scans at constant elevation. This scan strategy, combined with sky rotation, 
produced a W-like scan pattern in equatorial coordinates (Figure 1 of deB94), sparsely sampling 
a region 1 hour wide in right ascension and 15° wide in declination. deB94 remove a single 
offset from the entire data set prior to binning. The ARGO Hercules scan absolute calibration 
uncertainty is 5% (one standard deviation)^]. 

ARGO data is able to robustly constrain the normalization of the CMB spectra. However, 
ARGO alone is unable to strongly constrain the other cosmological parameters considered. Robust 



The value quoted in deB94, 15%, is a higher standard deviation calibration uncertainty. 
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constraints on these cosmological parameters will require a combined likelihood analysis of many 
CMB anisotropy data sets (including ARGO). In this paper we use the ARGO data set as a basis 
for developing a formalism that will eventually allow such a combined analysis. 

In §2 we summarize the computational techniques used in our analysis. See GRGS for a more 
detailed description. In §2 we also describe the greatly extended cosmological parameter space we 
consider here. Results are given and discussed in §3. Conclusions are presented in §4. 

2. Summary of Computation and Models Considered 

The reduced 2 mm Hercules data is shown in Figure 1. The zero-lag window function for 
these observations is shown in Figure 2 and the zero-lag window function parameters are listed in 
Table 1. 

Figure 2 also shows some of the model CMB anisotropy spectra used in the analysis. In earlier 
analyses (GRGS; Ganga et al. 1997b, 1998; Ratra et al. 1998), only a few (25) spectra were used. 
In this paper we use 798 spectra to cover a greatly extended range of the cosmological parameter 
space of the open and spatially-flat-A CDM models, and cover it with a higher resolution grid 
than previously used. This allows for an explicit construction of the full likelihood function as 
a function of all the cosmological parameters considered, and for making the various likelihood 
functions derived by marginalizing over parameters. This had not previously been possible. 

We focus here on a spatially open CDM model and a spatially flat CDM model with a 
cosmological constant A. These low density models are consistent with most current observational 
constraints. For discussions see Coles et al. (1998), Bartelmann et al. (1998), Jenkins et al. 
(1998), Park et al. (1998), Merchan et al. (1998), Cole et al. (1998), Cavaliere, Menci, k Tozzi 
(1998), and Somerville & Primack (1998). 

The models have gaussian, adiabatic primordial energy-density power spectra. The flat-A 
model CMB anisotropy computations use a scale-invariant energy-density perturbation power 
spectrum (Harrison 1970; Peebles & Yu 1970; Zel'dovich 1972), as predicted in the simplest 
spatially-flat inflation models (Guth 1981; Kazanas 1980; Sato 1981a, b). The open model 
computations use the energy-density power spectrum (Ratra & Peebles 1994, 1995; Bucher, 
Goldhaber, & Turok 1995; Yamamoto, Sasaki, & Tanaka 1995) predicted in the simplest 
open-bubble inflation models (Gott 1982; Guth & Weinberg 1983). The computation of the CMB 
anisotropy spectra is described in Stompor (1994) and Sugiyama (1995). 

The spectra are parameterized by their quadrupole- moment amplitude Qrms-PS, the 
clustered-mass density parameter Qq, the baryonic-mass density parameter Qb, and the age of 
the universe to- We have evaluated the spectra for a range of f^o spanning the interval 0.1 to 1 in 
steps of 0.1, for a range of f2#/i 2 [the Hubble parameter h = Ho/(WOh km s _1 Mpc -1 )] spanning 
the interval 0.005 to 0.029 in steps of 0.004, and for a range of to spanning the interval 10 to 20 
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Gyr in steps of 2 Gyr. Current observational estimates are more consistent with f2o ~ 0.4 (e.g., 
Cole et al. 1997; Eke et al. 1998), with t ~ 12 Gyr (e.g., Reid 1997; Gratton et al. 1997), and 
disfavor Q B h 2 larger than ~ 0.02 (e.g., Fukugita, Hogan, & Peebles 1998; Buries & Tytler 1998). 

While it is of interest to also consider other cosmological parameters, such as tilt or gravity 
wave fraction, to make the problem tractable we have focussed on the four parameters mentioned 
above. We emphasize however that the results of the analysis are model dependent. For instance, 
tilted flat- A models at fixed f2#/i 2 would likely lead to a different constraint on fio than that 
derived below in the scale- invariant flat- A model with varying fifi/i 2 . 

GRGS describe the computation of the likelihood function. Since deB94 remove an offset 
from the data, we assume a uniform prior in the amplitude of the offset removed and marginalize 
over this amplitude when computing the likelihood function (Bond et al. 1991, Bunn et al. 1994; 
GRGS; Church et al. 1997; Ratra et al. 1998). This must be done since the removal of the offset 
also removes an undetermined amount of the CMB anisotropy signal from the data. Beamwidth 
and calibration uncertainties are accounted for as described in GRGS. 

The open and flat-A model likelihoods are a function of four parameters: Qrms-PSj ^sh 2 , 
and to. We also compute marginalized likelihood functions by integrating over one or more of 
these parameters after assuming a uniform prior in the relevant parameters. The prior is set to 
zero outside the ranges considered for the parameters. 

To derive central values and limits from the likelihood functions we assume a uniform prior 
in the relevant parameter, so the corresponding posterior probability density distribution function 
vanishes outside the chosen parameter range and is equal to the likelihood function inside this 
range. The deduced central value of the parameter is taken to be the value at which the posterior 
probability density peaks. The limits we quote are based on the highest posterior density (HPD) 
prescription. They are determined by integrating the posterior probability density starting at the 
peak and minimizing the difference between the upper and lower limits. The ±1 a and ±2 a HPD 
limits encompass 68.3% and 95.5% of the area, respectively. See GRGS for further details. Of 
course, the quoted limits depend on the prior range considered for the parameter, if the likelihood 
function is not sharply peaked within the range considered. This is the case for a number of the 
likelihood functions obtained below. 



3. Results and Discussion 

For the flat bandpower spectrum the ARGO likelihood peaks at bandtemperature 5Ti = 33 [iK, 
with ala range of 28 < 5T t < 38 fiK and likelihood ratio = 1 x 10 21 . For the fiducial CDM 
model spectrum 8T\ = 30 fiK, with 1 a range 26 /iK < 57} < 35 fiK and likelihood ratio 5 x 10 21 . 
These numerical values account for the marginalization over the offset removed from the data, 
as well as the beamwidth and calibration uncertainties. For fiducial CDM deB94 find that the 
likelihood peaks at 57} = 39 fiK, with 1 a range 34 fiK < 57} < 45 and likelihood ratio = 
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5 x 10 23 . The deB94 1 a range accounts for calibration uncertainty by adding it in quadrature 
to the error bars derived ignoring it, but ignores offset removal and beamwidth uncertainty. The 
difference between the numerical values derived in deB94 and here is mostly because they have 
not marginalized over the amplitude of the offset removed while we have. 

For the flat bandpower spectrum the ARGO average 1 a STi error bar is ~ 15%: ARGO data 
results in a very significant detection of CMB anisotropy, even after accounting for beamwidth 
and calibration uncertainties. For comparison, the corresponding DMR error bar is ~ 10 — 12% 
(depending on model, Gorski et al. 1998). The MAX 4+5 data set also results in a small error bar 
~ 14% (Ganga et al. 1998). 

We note that <5TJ values estimated using the flat bandpower and fiducial CDM spectra are 
not identical. The variation in the deduced 5Ti values from model to model is an indication of the 
accuracy of the flat bandpower approximation to the real spectrum over the range of the window 
function of the experiment. The variation from model to model found here is comparable to that 
found from the SP94 and MAX 4+5 data (GRGS; Ganga et al. 1998) but is smaller than that 
found from the SuZIE and White Dish data (Ganga et al. 1997b; Ratra et al. 1998). 

For both the open and flat-A models, the four-dimensional posterior probability density 
distribution function L(Q rms _ps, £lo, ^_b/i 2 , to) is nicely peaked in the Qrms-PS direction but quite 
flat in the other three directions. The dotted lines in Figure 3 illustrate this flatness in the (f2o> *o) 
subspace of the flat-A model. Note that it is possible to distinguish between regions of parameter 
space at only slightly better than 0.25 a confidence. The irregular solid lines in Figure 5 are the 
2 a contours of the four-dimensional posterior distribution projected on to the (Oo, Qmis-Ps) 
subspace. They clearly show the steepness and peak (solid circles) in the Q rm s-PS direction. 

Marginalizing over Qmis-PS results in a three-dimensional posterior distribution 
L(Qq, Og/i 2 , to) which is steeper. The dashed lines in Figure 3 are the contours of this 
function. Note that it is now possible to distinguish between regions of parameter space at better 
than 1 a confidence. Marginalizing over one more parameter (in addition to Q rm s-Ps) results in 
two-dimensional posterior probability distribution functions which are peaked (albeit, in most 
cases, at an edge of the parameter range considered, see the solid circle in Figure 3). It is now 
possible to distinguish between parts of parameter space at better than 3 a confidence — see the 
solid lines in Figure 3. 

As discussed below, caution must be exercised when interpreting the discriminative power of 
these formal limits, since they depend sensitively on the values of the parameters beyond which 
the uniform prior has been set to zero. 

Figure 3 also illustrates a point noted earlier in the analysis of the MAX 4+5 data set (Ganga 
et al. 1998): conclusions about the most favored model drawn from the full four-dimensional 
posterior distribution tend to differ from those deduced from the three-dimensional posterior 
distribution derived by marginalizing the four-dimensional one over Q rm s-PS- For example, 
from a similar plot of the (SIq, ^b/j 2 ) subspace of the open model (not shown), one finds that 
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the four-dimensional distribution favors a large value of £Iq ~ 1, while the three-dimensional 
distribution favors a small f^o ~ 0.25. As discussed in Ganga et al. (1998), this is a consequence 
of the asymmetry of the posterior distribution. This is not a significant issue since at the 2 a 
confidence level the four-dimensional posterior distribution is flat in the Oo direction and so it is 
not statistically meaningful to discuss how it varies with Qq. 

As mentioned above, the two-dimensional posterior distributions (derived by marginalizing 
the four-dimensional distribution over two parameters at a time) allows one to distinguish between 
regions of parameter space at a higher level of confidence. Figure 4 illustrates this for the three 
cosmological parameters (Qo, f2_e/i 2 , io)> f° r both the open and flat-A models. Certain parts of 
parameter space can now be formally ruled out at better than 2 a significance. For example, for 
the open model a region in parameter space centered near Q ~ 0.7, S7#/i 2 ~ 0.005, and t ~ 10 
Gyr is ruled out at 3 a. Again, caution must be exercised when interpreting the discriminative 
power of such formal limits. 

Figure 5 shows the contours of the two-dimensional posterior distribution for Q r ms-PS an d 
fio, derived by marginalizing the four-dimensional distribution over fig/i 2 and to- These are shown 
for the ARGO, DMR, SP94, MAX 4+5, White Dish, and SuZIE data sets, for both the open and 
flat-A models. For all but the DMR data we also show 2 a confidence contours determined from 
projecting the four-dimensional posterior distribution on to this (Q r ms-PS> ^o) subspace. Clearly, 
at 2 <7, constraints on these parameters derived from the ARGO data are mostly consistent with 
those derived from the other data sets. However, the DMR and ARGO data are more consistent 
for an open model than for the flat-A case, panels o) and b) of Figure 5. In fact, there is very 
little overlap between the 2 a ranges (derived from the two-dimensional posterior distributions) of 
the ARGO and DMR normalizations for the flat-A models. There is, however, significant overlap 
between the 2 a ranges derived by projecting the ARGO and DMR four-dimensional posterior 
distributions^]. The SP94 results shown in panels c) and d) are those derived from the full Ka+Q 
data set. Sample variance and noise considerations indicate that the SP94 Ka band data is more 
consistent with what is expected for CMB anisotropy data (GRGS). At fixed Oo the Ka band 
data results in a lower deduced Q rm s-PS (compared to that from the Ka+Q data, GRGS), so the 
SP94 Ka band results are more consistent with the ARGO results. The MAX 4+5 results shown 
in panels e) and /) of Figure 5 are those derived from the MAX 4 ID and SH and MAX 5 HR, 
MP, and PH data sets. Sample variance and noise considerations indicate that the MAX 4 ID and 
MAX 5 HR data are more consistent with what is expected for CMB anisotropy data (Ganga et 
al. 1998). This subset of the MAX 4+5 data has a lower Q rms _ps (Ganga et al. 1998) and is thus 
more consistent with the ARGO data. 

The Python (Piatt et al. 1997) and Saskatoon (Netterfield et al. 1997) experiments are 



10 Since the DMR four-dimensional posterior distribution is only very weakly dependent on fish 2 and to, the 
confidence contours derived from this distribution projected on to the (Qrms-ps, Slo) plane are close to those derived 
from this distribution marginalized over Qsh 2 and to- 
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also sensitive to angular scales probed by ARGO. The Python large-chop result and some of the 
relevant Saskatoon n-point chop results are larger than the ARGO result. On the other hand, 
some of the Saskatoon n-point chop results are consistent with the ARGO result. 

Figure 6 shows the one-dimensional posterior distribution functions for S7o, &Bh 2 , to, and 
Qrms-PS) derived by marginalizing the four-dimensional posterior distribution over the other 
three parameters. From these one-dimensional posterior distributions, ARGO data favors an 
open (flat-A) model with either ft = 0.23 (0.10), or Q B h 2 = 0.029 (0.020), or t = 20 (12) Gyr, 
amongst the models considered. We emphasize that each of these are derived from one-dimensional 
posterior distributions and thus can not be simultaneously imposed. Also shown in Figure 6 are 
the limits derived from the one-dimensional distributions and the (projected) four-dimensional 
distribution. At 2 a confidence the ARGO data only formally rules out small regions of parameter 
space. Specifically, from the one-dimensional posterior distributions, the ARGO data requires f2o 
< 0.69 or > 0.76 (U < 0.91), or Q B h 2 > 0.006 (n B h 2 > 0.006), or t > 10 Gyr (t < 20 Gyr) 
for the open (flat-A) model at 2 a. Since some of the one-dimensional posterior distributions 
peak at the edge of the parameter range considered some of these limits must be considered to be 
formal. Less controversially, it is clear that, for both open and flat-A models, ARGO data favors a 
low-density universe, although Qq = 1 is not strongly ruled out. Similarly, ARGO data also mildly 
favors a relatively large Vt B h 2 or a young (old) flat-A (open) model. As discussed above, ARGO 
results in fairly tight constraints on Q rms -ps (panels g) and h) of Figure 6), and these are more 
consistent with the DMR results for the open model than for the flat-A case. 

Care is needed when interpreting the discriminative power of such formal limits. Consider 
a posterior density function which is a gaussian and nicely peaked inside the parameter range 
considered. The 1 and 2 a HPD limits for such a gaussian correspond to a value of the posterior 
distribution relative to that at the peak of 0.61 and 0.14 respectively. For the open model posterior 
distributions shown in Figure 6, the 1 and 2 a HPD limits correspond to values of the posterior 
distribution relative to that at the peak of 0.63 & 0.56 (Q ), 0.95 & 0.91 (tt B h 2 ), 0.93 & 0.86 
(to), and 0.36 & 0.13 (Qrms-ps)- It is hence probably fairer to conclude that the formal ARGO 
statistical limits on £l B h 2 and to, and the 2 a limits on S7o, should be taken much less seriously 
than those on Qrms-PS, and the 1 a limits on Qq. 

4. Conclusion 

In our likelihood analyses of the ARGO Hercules scan CMB anisotropy data we have explicitly 
accounted for beamwidth and calibration uncertainties and have marginalized over the amplitude 
of the offset removed from the data. As a consequence the results derived here differ from those 
derived earlier. The ARGO results are mostly consistent with those derived from the DMR, SP94, 
MAX 4+5, White Dish and SuZIE data sets. We emphasize again that the ARGO results, as well 
as those derived for the other data sets, depend on the models considered. 
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While the ARGO data does significantly constrain Qi-ms-PS ( an d mildly constrain Qq), robust 
constraints on these and other cosmological parameters from the CMB anisotropy must await a 
models-based combined likelihood analysis of many different data sets. 
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work was partially carried out at the Infrared Processing and Analysis Center and the Jet 
Propulsion Laboratory of the California Institute of Technology, under a contract with the 
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EPS-9550487 with matching support from the state of Kansas and from a K*STAR First award. 
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grant 2P03D00813. 



Table 1: Numerical Values for the Zero-Lag Window Function Parameters* 



l e -°- 5 l e -0.5 y/I(Wi) 

60 97.6 109 168 0.551 



a The value of I where Wi is largest, i m , the two values of I where Wi _ s = e 0,5 Wi m , i e -o.s, the effective multipole, 
lc = I(lWi)/I(Wi), and /(Wi) = XX 2 (Z + 0.5)Wi/{J(J + 1)}. 
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Figure Captions 

Fig. 1.- Measured 2 mm thermodynamic temperature differences (with ±1-<t error bars) on the 
sky as a function of scan position in the direction of Hercules. Note that the 63 points in the scan 
are at varying right ascension and declination. 

Fig. 2.- CMB anisotropy multipole moments 1(1 + l)Ci/(2ir) x 10 10 (solid lines, scale on left 
axis, note that these are fractional anisotropy moments and thus dimensionless) as a function of 
multipole I, for selected models normalized to the DMR maps (Gorski et al. 1998; Stompor 1997). 
Panels a) — c) show selected flat- A models. The heavy lines are the $7o = 0.1, Q B h 2 = 0.021, and 
to = 12 Gyr case, which is close to where the likelihoods (marginalized over all but one parameter 
at a time) are at a maximum. Panel a) shows five Q B h 2 = 0.021, to = 12 Gyr models with S7o = 

0. 1, 0.3, 0.5, 0.7, and 0.9 in descending order at the I ~ 200 peaks. Panel b) shows seven ^o = 0.1, 
t = 12 Gyr models with Q B h 2 = 0.029, 0.025, 0.021, 0.017, 0.013, 0.009, and 0.005 in descending 
order at the I ~ 200 peaks. Panel c) shows six Qq = 0.1, Q, B h? = 0.021 models with to = 20, 18, 
16, 14, 12, and 10 Gyr in descending order at the I ~ 200 peaks. Panels d) — f ) show selected open 
models. The heavy lines are the Qq = 0.2, Q B h 2 = 0.029, and to = 20 Gyr case, which is close 
to where the likelihoods (marginalized over all but one parameter at a time) are at a maximum. 
Panel d) shows five Q B h 2 = 0.029, t = 20 Gyr models with Q = 1, 0.8, 0.6, 0.4, and 0.2 from 
left to right at the peaks (the peak of the fio = 0.2 model is off scale). Panel e) shows seven $7o 
= 0.2, t = 20 Gyr models with Q B h 2 = 0.029, 0.025, 0.021, 0.017, 0.013, 0.009, and 0.005 in 
descending order at / ~ 600. Panel /) shows six Uq = 0.2, Q B h 2 = 0.029 models with to = 20, 18, 
16, 14, 12, and 10 Gyr in descending order at I ~ 600. Also shown is the ARGO 2 mm zero-lag 
window function W\ (dotted lines, scale on right axis). See Table 1 for ^-parameter values. The 
ARGO data mainly constrains the area under the product of the model spectrum and the window 
function on this plot. This is clearly larger for the DMR-normalized flat- A models than for the 
DMR-normalized open models. ARGO thus favors a lower relative normalization between the 
flat-A and open model than does the DMR. 

Fig. 3.- Maxima and confidence contours of various posterior probability density distribution 
functions for the (f2o,to) subspace of the flat-A model. Dotted lines (solid squares) show the 
contours (maxima) of the four-dimensional (Q r ms-PS> ^Bh 2 , to) posterior distribution; contours 
of 0.1 and 0.25 a confidence are shown. Dashed lines (solid triangles) show the contours (maxima) 
of the three-dimensional (Q,Q,Q, B h? , to) posterior distribution (derived by marginalizing the 
four-dimensional one over Q rm s-ps)- Contours of 0.1, 0.25, 0.5, and 1 a confidence are shown. 
Solid lines (solid circles) show the contours (maxima) of the two-dimensional posterior distribution 
(derived by marginalizing the four-dimensional one over the other two parameters). Contours of 

1, 2, and 3 a confidence are shown (the 3 a contour is not labelled). 

Fig. 4.- Confidence contours and maxima of the two-dimensional posterior probability density 
distribution functions, as a function of the two parameters on the axes of each panel (derived by 
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marginalizing the four-dimensional posterior distribution over the other two parameters). Dashed 
lines (open circles) show the contours (maxima) of the open case and solid lines (solid circles) 
show those of the flat-A model. Contours of 0.25, 0.5, 1, 2, and 3 a confidence are shown (3 a 
contours are not labelled). 

Fig. 5.- Confidence contours and maxima of the two-dimensional (Qrms-PSi ^o) posterior 
probability density distribution function. Panels a), c), e), and g) are for the flat-A model and 
panels b), d), /), and h) are for the open model. Note the different scale on the vertical (Qrms-ps) 
axis in each pair of panels. Shaded regions show the ARGO results, with denser shading for the 
1 a confidence region and less-dense shading for the 2 a region. Irregular solid lines show the 2 
a confidence contours derived by projecting the four-dimensional ARGO posterior distribution in 
to this plane. [I.e., for each set of values of (<2rms-PS> ^o) we check if there is any choice of Sl^/i 2 
and to such that the point (Q r ms-PS; ^o> &Bh 2 , to) is within the 2 a region of the four-dimensional 
ARGO posterior distribution. The projected 2 a limits enclose those values of (Qrms-PSi ^o) f° r 
which such a point exists.] In panels a) — f) hatched areas show the two-dimensional posterior 
probability density distribution function confidence regions for the DMR data (panels a) and b), 
Gorski et al. 1998; Stompor 1997), the SP94 Ka+Q data (panels c) and d), GRGS), and the MAX 
4+5 data (panels e) and /), Ganga et al. 1998). Heavy dashed confidence contours bounding 
these regions are labelled (except for the DMR cases); denser hatching corresponds to the 1 a 
confidence region and less-dense hatching to the 2 a confidence region. In panels c) — /) unlabeled 
irregular light dashed lines show the 2 a confidence contours derived by projecting the SP94 and 
MAX 4+5 four-dimensional posterior distributions in to this plane. Panels g) and h) show the 
SuZIE 2 a upper limit (the hatched region bounded by the labelled heavy dashed line, Ganga et 
al. 1997b) and the White Dish 2 a upper limit (labelled heavy dotted line in panel h) and not 
shown in panel g) since it is off scale, Ratra et al. 1998). The unlabeled irregular light dashed 
and dotted lines in panel h) are the corresponding SuZIE and White Dish 2 a confidence upper 
limits derived by projecting the four-dimensional posterior distribution in to this plane (these 
limits are not shown in panel g) since they are off scale). Solid circles show the maxima of the 
ARGO two-dimensional posterior distribution and open circles show those of the other data sets 
(not shown for DMR, SuZIE, and White Dish). The DMR results are a composite of those from 
analyses of the two extreme data sets: i) galactic frame with quadrupole included and correcting 
for faint high-latitude galactic emission; and ii) ecliptic frame with quadrupole excluded and no 
other galactic emission correction (Gorski et al. 1998). 

Fig. 6.- One-dimensional posterior probability density distribution functions for S7o, &Bh 2 , to, 
and Qrms-PS (derived by marginalizing the four-dimensional one over the other three parameters) 
in the open and flat-A models. These have been renormalized to unity at the peaks. Dotted 
vertical lines show the confidence limits derived from these one-dimensional posterior distributions 
and solid vertical lines in panels g) and h) show the ±1 and ±2 a confidence limits derived by 
projecting the four-dimensional ARGO posterior distribution. Note that, as discussed at the end of 
§3, some of these formal limits have very little discriminative power. The 2 a DMR (marginalized 
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and projected) confidence limits in panels g) and h) are a composite of those from the two extreme 
DMR data sets (see caption of Figure 5). 
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